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ABSTRACT 

The early reionisation of the Universe inferred from the WMAP polarisation results, 
if confirmed, poses a problem for the hypothesis that scale-invariant adiabatic density 
fluctuations account for large-scale structure and galaxy formation. One can only gen- 
erate the required amount of early star formation if extreme assumptions are made 
about the efficiency and nature of early reionisation. We develop an alternative hy- 
pothesis that invokes an additional component of a non-scale-free isocurvature power 
spectrum together with the scale-free adiabatic power spectrum for inflation-motivated 
primordial density fluctuations. Such a component is constrained by the Lyman alpha 
forest observations, can account for the small-scale power required by spectroscopic 
gravitational lcnsing, and yields a source of early star formation that can reionise the 
universe at z ~ 20 yet becomes an inefficient source of ionizing photons by z ~ 10, 
thereby allowing the conventional adiabatic fluctuation component to reproduce the 
late thermal history of the intergalactic medium. 
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1 INTRODUCTION 

There have been major recent advances in cosmology with 
impact on galaxy formation theory. These include the detec- 
tion of the temperature-polarisation cross power spectrum 
for the CMB by WMAP (Kogut et al. 2003), and measure- 
ments both of the CMB temperature power spectrum and 
the underlying matter power spectrum with unprecedented 
accuracy, utilizing the WMAP, 2DF and quasar Lya absorp- 
tion line data sets (Bennett et al. 2003; Spergel et al. 2003). 
Two results that have received considerable attention are 
the optical depth of the universe r = 0.17 ± 0.03, which re- 
quires that the epoch of reionisation occurs at z= 15-20 from 
WMAP (Kogut et al. 2003; Spergel et al. 2003), and the 
rolling spectral index dn/dlnk = —0.03 ± 0.01 for approxi- 
mately scale-invariant density fluctuations for a combination 
of WMAP, 2DF and Lya data (Spergel et al. 2003). 

There is some tension between these results: if both 
are correct, it is difficult to understand how recombina- 
tion occurred so early without some modification of the 
canonical model of primordial, nearly scale-invariant Gaus- 
sian adiabatic density fluctuations (Ciardi, Ferrara & White 
2003; Fukugita & Kawasaki 2003; Haiman & Holder 2003; 
Somerville et al. 2003). In fact, a new Lya absorption data 
set from the SDSS quasars has independently found evidence 



for a rolling spectral index (Seljak 2003), although an inde- 
pendent analysis of the same data does not reproduce suf- 
ficiently small error bars to confirm this result (Abazajian, 
& Dodelson 2003). 

The Lya lines measure power in the underlying mat- 
ter power spectrum on a comoving scale of around 1 Mpc. 
The results are however subject to bias, since one has to be 
confident that the gas is relatively unperturbed by feedback, 
such as is seen in the vicinity of Lyman break galaxies to 
Mpc distances. Hence it is of particular interest to consider 
another measure of the power spectrum on even smaller co- 
moving scales, 10 -2 to 0.1 Mpc. This comes from spectro- 
scopic gravitational lensing of quasar emission line region on 
several scales, the magnification ratios requiring and con- 
straining substructure in the massive lensing halos (Metcalf 
et al. 2003). One needs substantial power in objects of 10 6 to 
10 9 M©, amounting to between 4 and 7 percent of the galaxy 
surface density, and this cannot easily be accommodated 
in the usual CDM models with standard elliptical isother- 
mal lens mass profiles. Previous estimates of halo substruc- 
ture from gravitational lensing using simple lens models are 
highly uncertain (Dalai & Kochanek 2002). Moreover the 
numerical simulations of halos suggest that, for nearly scale- 
invariant initial conditions, small-scale power may largely be 
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erased in the inner halos. One is in a dangerous regime of 
the spectrum where n e s ~ —3, and tidal destruction is ef- 
fective. Indeed it is not clear whether the simulations have 
enough resolution to adequately address the question of the 
survival of small-scale substructure. 

We propose an alternative prescription for small-scale 
power that satisfies all observed constraints and unambigu- 
ously predicts the survival of small-scale power, with some 
inflationary motivation. We postulate that as is common to 
multifield inflationary models, both isocurvature and adia- 
batic fluctuation modes are generated. A recent model mo- 
tivated by inflation is a so-called curvaton model (Lyth, & 
Wands 2002; Moroi & Takahashi 2001; Enqvist & Sloth 

2002) in which an additional scalar field besides the inflaton, 
- the curvaton-, produces fluctuations during the reheating 
epoch. If the isocurvature fluctuations were generated by 
inflaton and curvaton generated adiabatic mode, there may 
be a possibility of having a non-scale-free isocurvature power 
spectrum together with the n = 1 scale- free adiabatic power 
spectrum (Moroi 2003) . As well as this curvaton hypoth- 
esis, a sub-dominant contribution of cosmic strings induced 
by brane inflation in superstring theory(Jones, Stoica & Tye 

2003) may provide additional small-scale power. 

The isocurvature contribution is described by two pa- 
rameters: the amplitude normalisation and the spectral in- 
dex, which we take to be freely assignable, but chosen to 
give more small-scale power than the rolling or nearly scale 
invariant index adiabatic fluctuations measured on larger 
scales. We use Lya forest data to constrain these parameters. 
The amount of excess small-scale power can be tuned by 
adjusting the spectral index within the allowed constraints. 
We normalise at 1 Mpc, the central point of the Lya probes. 
Small-scale power survives because the fluctuations become 
nonlinear earlier than in the pure adiabatic case, due to the 
isocurvature component boost. The first nonlinear fluctua- 
tions form earlier and hence lead to denser substructures 
that are resistant to tidal disruption within massive halos. 

The two-component model has two advantages. It re- 
sults in early star formation, regardless of the spectral in- 
dex measured for the adiabatic component on large scales. 
Hence early reionisation can be achieved. It also preserves 
small-scale power as hierarchical clustering develops, in the 
form of dense 10 6 Mq clumps in massive dark halos. This 
helps to explain quadruple quasar lensing flux ratios as well 
as bending of radio minijets. 

In the remainder of this paper, we give constraints from 
the Lya forest on the isocurvature component and we discuss 
the observational implications. We give several applications: 
in addition to the implications for the epoch of reionisation 
and halo microlensing, we discuss the possible implications 
for patchy reionisation and SZ signatures of very early star 
formation via baryon trapping in dense early substructures, 
and the clustering of early forming substructures and impli- 
cations for the formation of the first stars. 

2 MODELS 

From WMAP results, it is known that the nature of 
fluctuations is consistent with the adiabatic mode and 



the contribution of isocurvature perturbations cannot be 
dominant on the relevant measurement scales which are 
k < 0.05Mpc _1 (Spergel et al. 2003; Bennett et al. 2003). 

On smaller scales, k ~ IMpc -1 , the amplitude of isocur- 
vature fluctuations is constrained by Lya forest (Croft et 
al. f998; Nusser & Haehnelt 1999). The joint analysis of 
WMAP, 2dF and Lya shows the power spectrum obtained 
by Lya is significantly lower than the one extrapolated from 
a single power law A CDM model consistent with WMAP 
data alone (Spergel et al. 2003). 

Several new analyses of Lya power spectrum based 
on the Sloan Digital Sky Survey data are becoming avail- 
able (Seljak 2003; Abazajian, & Dodelson 2003). Although 
they employed the same data sets, each group has obtained 
different values. Among them, Seljak (2003) provides the 
lowest amplitude of the power spectrum and the smallest 
error bars. Therefore here we take his value as a reference. 
We set the normalisation of the isocurvature fluctuations to 
be 10% of Seljak's power spectrum. We should notice that 
this is merely an upper limit on the possible isocurvature 
contributions. 

On the other hand, Abazajian & Dodelson (2003) give 
a relatively high central value with larger error bars. The 
central value of their power spectrum indeed exceeds the 
extrapolation from the WMAP power law A CDM power 
spectrum. Here we have one extreme model assuming the 
difference between the Abazajian & Dodelson power spec- 
trum and the WMAP power law spectrum is due to the 
existence of the isocurvature contributions. Then we imme- 
diately obtain the amplitude and the power law index (which 
is n = —1.7) of the isocurvature power spectrum. The cor- 
responding power spectra are shown in Figure 1. 



3 REIONISATION 

We estimate the number of photons emitted from Pop III 
stars, which is a crucial element for the reionisation process, 
following Somerville & Livio (2003) and Somerville, Bullock 
& Livio (2003). 

First, employing the Press- Schechter prescription, we 
calculate the fraction of the total mass in collapsed halos Fh 
with masses greater than M cr it and lower than M V i r . Here 
we adopt M cirt = 1 x W fi h~ 1 M e and M vir = M(T vir = 
10 4 K). Objects whose virial temperature are larger than 
10 4 K can cool via atomic hydrogen and we assume them 
as Pop II stars. 

From Fh, we can calculate the global star formation rate 
density as 

dFh 

P* = e,p B —{M vir > M > M cri t) , (1) 

where e* is the star formation efficiency which we take 0.002 
for Pop III stars with 200M Q and 0.001 with 100M©, and 
pB is the comoving background baryon density (Yoshida et 
al. 2003). 

Let us assume Pop III stars produce dN^/dt = 1.6 x 10 48 
photons s _1 Mg 1 for a lifetime At — 3 x 10 years (Bromm 
et al. 2001). Therefore we can write dN 7 /dt = A 7i0 6(t), 
where A 7 ,o = 1.6 x 10 48 photons s _1 Mq 1 and O(t) is a step 
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Figure 1. The matter power spectra. The red line is the WMAP 
best fitted power law A CDM model (PLACDM) (Spergel et al. 
2003), the black hatched region is the running spectral index 
(RSI) model with errors fitted by WMAP, 2dF and Lyo (Spergel 
et al. 2003), and the blue and green hatched regions are the RSI 
model with the SDSS Lyo power spectrum analyzed by Seljak 
(2003), and Abazajian & Dodelson (2003), respectively. Purple 
lines, which are labeled as LFMIN (Lya Forest Minimum), are 
isocurvature power spectra with n = —2.5, —2, —1.5, and —1 nor- 
malized to 10% of Seljak's Lyo analysis. The pink line labeled as 
LFMAX (Lyct Forest Maximum) is the isocurvature power spec- 
trum by assuming the excess of the central value of Abazajian & 
Dodelson (2003) is caused by isocurvature contributions. 

function as Q(t) = 1 for t < At and @(t) = for t > At. 
Using this expression, the total production rate of ionizing 
photons per cubic Mpc at time t' becomes 

ft' jp 
= e tPB N lfi J Q( t ' -t)-^(t)dt 

f*' dF h , , , 

Jt'-At dt 

= e, PB N lfi (F h {t')-F h {t' -At)). (2) 

Then we can obtain the cumulative number of photons per 
H atom as 

cumul pt a , 

Lh fc,)=£™E / ^(t')dt' 

n H pB J dt 

rto 

=/im p e,iV 7 ,o / (F h (t') - F h (t' - At)) dt', (3) 
Jo 

where nu is the Hydrogen number density and m p is 



Figure 2. Cumulative photons emitted from Pop III stars. Mod- 
els are same as Figure 1 while we take the central value for the 
running spectral index model (RSI) and Abazajian, & Dodelson 
SDSS analysis (A-D). The difference between power law ACDM 
(PLACDM) and A-D is very small. 



the proton mass. Here we may employ the approxima- 
tion F h (t') - Fh(t' - At) ~ dF h /dt'(t')At, and we obtain 
n cumui / nH (t ) ~ iim v e*N^ fi F h (t )At. We checked this ap- 
proximation works almost perfectly well for z < 40 when 
the Hubble time is longer than At. 

For calculating the number of photons emitted from Pop 
II stars, we replace M v i r > M > M CT u of equation (1) with 
M > M vir and set e* = 0.1 and dN 1 /dt = 8.9x 10 46 photons 
s^Mg 1 for a lifetime At = 3 x 10 6 years (Somerville & Livio 
2003). 



4 RESULTS 

In Figure 2, the cumulative number of photons emitted from 
Pop III stars per hydrogen atom as a function of redshift is 
shown for the power law ACDM model with WMAP pa- 
rameters, the running spectrum index model, isocurvature 
models with the power law index n = —2.5, —2, —1, —1.5, — 1 
whose amplitudes are set as 10% of Seljak's analysis of Lya 
forest, the isocurvature model with n = —1.7 which explains 
the excess of the power spectrum obtained by the analysis 
of Lycv clouds by Abazajian & Dodelson against power law 
ACDM and the model fitted with Abazajian & Dodelson. 
Here, we take e* = 0.002. However we can simply renor- 
malise the value of y-axis for different values of e* . 

In Figure 3, the cumulative number of photons emitted 
from Pop II stars per hydrogen atom as a function of redshift 
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Figure 3. Cumulative photons emitted from Pop III stars. Mod- 
els are same as Figure 2. 



is shown. It is clear from Figures 2 and 3 that contributions 
from Pop II stars on early reionisation are almost negligible. 

To get the ionisation fraction, we have to multiply 
fescfion/Cdump where f esc is the escape fraction of pho- 
tons from the galaxy, fi on is the number of ionisations per 
UV photon, and C c i ump is the clumping factor of IGM, re- 
spectively. It is known that about 5 to 20 photons per H 
atom are required to achieve a volume-weighted ionisation 
fraction of 99% (Haiman, Abel & Madau 2001; Sokasian et 
al. 2003; Sokasian et al. 2003). We take 20 as our reference. 
In Figure 2, this number is plotted as the horizontal line. If 
we assume e* = 0.001 instead 0.002, 20 photons per H atom 
correspond to 40 photons in this figure which we also plot- 
ted as a horizontal line. The crossing of each line with this 
horizontal line gives the epoch of reionisation. The vertical 
line in this figure is z = 17 which is the reionisation epoch 
for WMAP with instantaneous reionisation. 

In Table 1, we summarize the "reionisation" epoch of 
each model. These numbers correspond to e* = 0.002. If we 
employ e* = 0.001 and assume the cumulative number of 
photons needed for reionisation is between 10 and 40 per H 
atom, the corresponding redshifts at which this occurs are 
given in this Table. It is clear that the running spectrum 
index model cannot plausibly have early enough reionisa- 
tion and that power law ACDM is marginally consistent 
with WMAP results. These results are consistent with pre- 
vious works (Yoshida et al. 2003; Ciardi, Ferrara & White 
2003; Fukugita k, Kawasaki 2003; Haiman & Holder 2003; 
Somerville & Livio 2003) It should be noticed that these 
redshifts for models with isocurvature fluctuations besides 
the A-D model are upper limits since we choose the ampli- 



Table 1. Redshifts at which the cumulative number of photons 
becomes 10, 20 and 40. Redshifts when the cumulative number 
of photons from PopII becomes equal to the one from PopIII are 
also shown. 
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tudes of power spectra to be as large as possible without 
violating Lya constraints. 



5 DISCUSSION 

We can certainly have early enough reionisation if we intro- 
duce the isocurvature mode on small scales. One interesting 
aspect is the number of cumulative photons per H atom by 
PopIII stars asymptotes to a constant at a later epoch for 
some models. To illustrate this effect, we plot the number of 
cumulative photons for the n = — 1 model in Figure 4. Here 
we renormalise the power spectrum to have an appropriate 
reionisation epoch. We have the number of photons about 20 
at z = 17 for assuming 1% amplitude of Seljak's Lya power 
spectrum and 10 for 1%. We can clearly see the flattening of 
number of photons at z 7. This flattening may allow the 
small neutral H fraction which is observed in the spectra 
of the highest redshift SDSS QSOs, via the Gunn-Peterson 
effect. 



5.1 Halo microlensing 

The flux ratios of several quadruple-lensed quasars can only 
be interpreted if halo substructure is adding differentially 
to the lensing optical depth. Between 0.6 and 7 percent of 
the halo mass is required to be in structures of mass up to 
10 8 — 10 10 Mq, within a projected radius of 10 kpc of a mas- 
sive halo at z = 0.6 citedal,met. Most of the contribution 
to the optical depth comes from within the scale radius of 
the dark halo, since at larger radii the mean halo density 
decreases as r~ 3 . However the numerical simulations do not 
have the resolution to tell whether the halo substructure 
survives, for a canonical scale-invariant initial spectrum of 
fluctuations. Semi-analytical methods suggest that the sub- 
structure fraction is insensitive to tilt or roll, but possibly 
too low (Zentner and Bullock 2003) for the purely adiabatic 
model. 

The model advocated here can readily accommodate 
the needs of halo substructure lensing, as the early forming 
substructures are more numerous and denser, and so resis- 
tant to tidal disruption. 
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Figure 4. Cumulative photons emitted from Pop III stars with 
smaller normalisation for the n = — 1 isocurvature model. Only 
1% and 2% of Seljak's analysis of Lya forest are needed for 10 
and 20 photons par H atom at z = 17. 



5.2 The mass fraction in minihalos, Populations 
II and III at high z 

We may define minihalos to be dark matter clouds which are 
below the mass threshold for star formation. The relevant 
mass range for minihalos that can trap baryons requires tem- 
peratures above that of the CMB and masses above about 
1O 4 /i -1 M0. In contrast, cooling is only effective at masses 
above approximately 1O 6 M0. 

The abundance of minihalos is shown in Figure 5 as 
a function of redshift in a typical isocurvature/adiabatic 
model. Note that they are more numerous out to z ~ 40 
than the peak in the ACDM model, which occurs for mini- 
halos at z ~ 10. 

Also shown in Figure 5 are the mass fractions in Pop- 
ulation III and in Population II stars that form in dwarf 
galaxy halos. These are defined by the respective criteria 
that cooling by H2 and Lya cooling are the dominant dis- 
sipative mechanisms for concentrating the baryons and en- 
abling fragmentation to proceed. We base our criteria for 
formation of Population III and Population II stars in pri- 
mordial clouds on the formulation by Haiman and Holder 
(2003) in terms of Type II vs Type I halos. Their classifica- 
tion is based on the distinction between H2 and HI cooling: 
we simply take this definition to its logical conclusion, given 
that the consensus view is that molecular cooling results 
in very massive stars (Pop III) and atomic cooling allows 
fragmentation to the "normal" mass range (Abel, Bryan & 
Norman 2002; Bromm, Kudritzki & Loeb 2001; Omukai 
& Yoshi 2003). The corresponding mass ranges are defined 



Figure 5. The mass fraction of halos, within a given mass range, 
is shown as a function of redshift for the PLACDM model (solid 
lines) and the n = —1.5 isocurvature reionisation model (dashes 
lines). For each model 3 cases are presented: mass intervals defined 
by all masses with T v i r > 10 4 K (predominantly halos that cool 
via Lya and form popll stars); masses with T v i r < lO 4 -^ and 
above 10 6 h~ 'Mq (predominantly halos that cool via H2 and form 
popIII stars); and minihalos of mass above 10 4 h~ x Mq and less 
than 10 e h~ 'Mq (halos that are too low in mass to cool) but may 
contain residual trapped gas. 

by all masses with T vir < lO 4 -^ and above 1O 6 /i _1 M0 (pre- 
dominantly halos that cool via H2 and form popIII stars), 
and by all masses with T v t r > W 4 K (predominantly halos 
that cool via Lya and form Population II stars). 

We see that Population III stars are boosted by an or- 
der of magnitude at z ~ 20, although Population II star 
formation is not greatly affected by the isocurvature admix- 
ture relative to ACDM. This is because the mass fraction 
in the relatively massive clouds required in this latter case, 
typically in excess of ~ 10 9 /t -1 Mo, is strongly constrained 
by our model which incorporates the requirement that we 
cannot overly perturb the Lya forest. 



5.3 Baryon trapping and SZ fluctuations 

We consider the effects of baryon trapping in the isocur- 
vature perturbation-induced substructure. This will have 
the effect of enhancing the temperature and SZ fluctuations 
produced at reionisation relative to those predicted for the 
pure adiabatic case. Even if the baryons cannot cool, they 
are trapped at high redshift in dark matter minihalos of 
mass above about 1O 4 /i _1 M0. The baryon overdensity is 
~ (a v /v s ) 2 , where a v is the velocity dispersion in the dark 
matter minihalo and v s is the gas sound velocity. Trapping 
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occurs only if T > Tomb, and this happens in the more 
massive minihalos e.g. above 1O 4 /i -1 M0. Cooling via H2 fur- 
ther enhances the gas density in late-forming minihalos. It 
has been argued that Population III stars form in such gas, 
whereas once the H2 is photo-dissociated , H cooling pre- 
dominates via Lya excitations and stars of lower mass can 
form as fragmentation continues to higher densities (Omukai 
&Yoshi 2003). 

However star formation is by no means guaranteed. 
Minihalos may retain gas supported in a stable configuration 
by dark matter self-gravity (Umemura & Ikeuchi 1986; Ger- 
hard & Silk 1996). In the present case, such minihalos could 
be very abundant at z > 20, and may provide a unique 
window on the dark ages of the early universe via radio 
and NIR observations of a diffuse background of redshifted 
21cm and Lya emission. For example with 100 Lya photons 
per baryon one might see at 2/4 a 1 percent contribution 
to the diffuse extragalactic background, which amounts to 
vi v « 10nwm _2 sr _1 , but could however be spectrally con- 
centrated in a feature with width Avjv ~ 0.1 associated 
with the epoch of reionisation. 
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